Higgs Modulated Reheating of RG improved Inflation in light of the Planck data 
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Within the framework of RG improved inflationary cosmology motivated by asymptotically safe 
gravity, we study the dynamics of a scalar field which can be interpreted as the Higgs field. The 
background trajectories of this model can provide sufficient inflationary e-folds and a graceful exit to 
a radiation dominated phase. We study the possibility of generating primordial curvature perturba- 
tions through the Standard Model Higgs boson. This can be achieved under finely tuned parameter 
choices by making use of the modulated reheating mechanism. The primordial non-gaussianity is 
expected to be sizable in this model. Though tightly constrained by the newly released Planck CMB 
data, this model provides a potentially interesting connection between collider and early universe 
physics. 
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I. INTRODUCTION 

Inflation is a successful model for the early universe, 
and manages to address several challenges of Big Bang 
cosmology [T'JS] (see [3HS] for early works). One of in- 
flation's crowning achievements is the prediction of a 
nearly scale-invariant primordial power spectrum which 
was verified to high precision by Cosmic Microwave Back- 
ground (CMB) observations [THS]. The success of infla- 
tion, however, is based on a series of assumptions includ- 
ing the existence of a sufficiently long period of quasi- 
exponential expansion realized by an as yet unobserved 
slow-roll scalar field. An important question to be ad- 
dressed is whether such a scalar field with an appropriate 
potential indeed exists in nature. 

The Higgs boson is a scalar field predicted by the Stan- 
dard Model (SM) of particle physics. Recent results from 
the ATLAS and CMS collaborations at CERN's Large 
Hadron Collider (LHC) have established the existence 
of a new resonance at 125 GeV, which is compatible 
with the Higgs boson with Scr significance for the sig- 
nal [ini E] • As the only existing scalar particle in the 
SM, it is natural to speculate that the Higgs field may 
be used to play a dual role as the infiaton. However, 
the corresponding Higgs energy scale is much too low 
compared with the typical infiationary scale. In order 
to solve this problem, a non-minimal coupling was intro- 
duced between the Ricci scalar and the Higgs field [T^fTS] 
(see also an earlier attempt in il4 ), and the associated 
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asymptotic freedom was studied in [15]. It was soon re- 
alized that such a non-minimal coupling leads to a uni- 
tarity violation energy scale that must be larger than the 
inflationary scale; otherwise the effective field description 
would fail (for example, see [TBl - B^ for detailed discus- 
sions). One thus faces challenges when constructing a 
suitable inflation model based on the Higgs field. 

One of the most challenging tasks in physics is the 
construction of a consistent UV-complete theory of grav- 
ity. Weinberg made an intriguing proposal in the 1970's: 
that the effective description of quantum gravity may 
be non-perturbatively renormalizable via the notion of 
asymptotic safety (AS) [531 [23] . In such a scenario the 
flows of the renormalization group (RG) approach a fixed 
point in the ultraviolet (UV) limit, and a finite dimen- 
sional critical surface of trajectories evolves to this point 
at short distance scales [25l [26] . Such a fixed point can 
be found in the Einstein-Hilbert truncation [27] , and the 
AS scenario has been studied extensively in the literature 
[25H55] (for recent reviews see [M1I55]). 

Recently, two of the present authors (YFC and DAE) 
proposed a cosmological model based on the concept 
of asymptotically safe gravity, in which the Higgs bo- 
son plays an important role during inflation [3B]. While 
there is as yet no explicit proof of the AS behavior, with 
its theoretical promise it warrants to explore whether 
the gravity theories under such a scenario are experi- 
mentally testable. It has been argued that if there are 
no intermediate energy scales between the SM and AS 
scales, the mass of the Higgs boson is predicted to be 
ruH = 126 GeV, with only several GeV uncertainty [37]. 
Inspired by this idea, in our model we consider the grav- 
itational constant G and the cosmological constant A 
varying with the energy scale, and thus varying through- 
out the cosmological evolution. Since the gravitational 
constant and the cosmological constant vary along with 



the running of the cutoff scale, the stress-energy tensor 
acquires an extra contribution. When taking into ac- 
count the extra term, the model is found to correspond 
to a f{R) model [5S] (see also [321 HO] for generalized 
discussions). As a consequence, there are effectively two 
scalar degrees of freedom, one being the adiabatic mode 
and the other being an iso-curvature mode. We find that 
the corresponding perturbation theory leads to primor- 
dial power spectra for both the curvature perturbation 
and the entropy perturbation. When the cutoff scale 
runs lower than certain critical value, the universe exits 
inflation gracefully and the AS gravity approaches GR 
simultaneously. 



II. ASYMPTOTICALLY SAFE GRAVITY 

We begin with a brief review of the cosmological of 
AS gravity minimally coupled to ordinary matter La- 
grangian. Specifically, we consider AS gravity with 
Einstein-Hilbert truncation and corresponding beta func- 
tions including next-to-leading order corrections of grav- 
itational coupling. This type of gravitational theory can 
be reformulated as an f{R) theory non-minimally cou- 
pled to the matter Lagrangian. 



A. Asymptotic safety 



After inflation, there exist two possible mechanisms of 
converting entropy fluctuations into curvature perturba- 
tions. The first is to make use of the curvaton mech- 
anism [HHHj . However, in the model of AS inflation 
such a transfer is highly nonlinear and thus leads to se- 
vere fine-tuning problem on model parameters [36 (see 
also [13 HTj for Higgs-curvaton problems in GR). Re- 
cently, a second possibility of realizing such a transfer 
was suggested in Refs.[i51 H5] . where the Higgs boson 
could modulate the efSciency of reheating |50n52] . This 
mechanism requires that the decay rate of the inflaton 
to be a function of a second field with a much lighter 
mass. This allows the entropy fluctuations during infla- 
tion to transfer to curvature perturbations on a hyper- 
surface with the average Hubble expansion comparable 
to the inflaton decay rate. 



In the present paper we study the possibility of AS in- 
flation with modulated reheating through the Higgs bo- 
son. The paper is organized as follows. Section H gives 
a brief review on the AS scenario of gravity theory with 
Einstein truncation and its classical correspondence to 
f{R) gravity is introduced. In Section HI, we perform 
a conformal transformation from Jordan frame to Ein- 
stein frame and then study the inflationary background 
driven by the RG improved gravitational and cosmolog- 
ical constants. We show that this cosmological system 
is equivalent to a double field inflation model where the 
contribution of the Higgs boson is secondary. Section IV 
is devoted to the analysis of modulated reheating via the 
Higgs field during which the isocurvature fluctuation is 
converted into the adiabatic perturbation and thus forms 
the primordial power spectrum of curvature perturba- 
tion. The primordial non-gaussianity seeded by nonlin- 
ear fluctuations of the Higgs is computed as well. In Sec- 
tion V, we confront our model with the newly released 
Planck data and find that the parameter space is tightly 
constrained. Finally, we conclude with a discussion in 
Section VI. We will work with the reduced Planck mass, 
Mp = I/^/SttGn, where Gn is the gravitational constant 
in the IR limit, and adopt the mostly-plus metric sign 
convention (— , +, +, +). 



We start with a RG inspired effective gravitational La- 
grangian with Einstein-Hilbert truncation. 
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where p is the RG cutoff scale beyond which the UV 
modes are argued to be integrated out. This effective La- 
grangian automatically connects with ordinary Einstein 
gravity in the IR regime where the gravitational and cos- 
mological constants flow to some integral constants that 
are constrained by observations. In the UV limit, these 
"constants" flow to a UV fixed point according to their 
beta functions. Quantum corrections are therefore en- 
coded in the evolution of the coupling constants as func- 
tions of the cutoff scale, whose beta functions can be 
extracted from the RG equations. 

We define the dimensionless gravitational and cosmo- 
logical constants as follows, 
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Given the exact forms of RG equations, one can follow 
the flows of g and A along with the cutoff scale and obtain 
a fixed point in the UV limit. The AS scenario suggests 
that this UV fixed point is attractive. Note that the ex- 
plicit forms of beta-functions depend on the choice of the 
cutoff function and the relevant gauges. In Ref. [5F, it 
was observed that the UV fixed point often corresponds 
to a de Sitter solution, however neither the energy scale 
of the background nor the amplitude of quantum fluctu- 
ations provide a successful application to early universe 
inflationary cosmology [38 . 

If the RG-improved gravity theory is viable, its RG tra- 
jectory should connect smoothly with standard Einstein 
gravity in the IR limit so as to be consistent with astro- 
nomical and cosmological observations. We would there- 
fore like to study the RG improved gravity theory in the 
regime that is sufficiently close to GR while still retain- 
ing linearized quantum corrections to the beta-functions. 
To begin the analysis, we linearize the beta-functions for 
dimensionless coupling constants as follows: 
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which include next-to-leading order corrections to g. The 
coefBcients a and /3 are cutoff functions dependent. Un- 
der this parameterization, one can obtain approximate 
forms of the dimensionless couplings, which are given by 
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where Gat and A/jj are Newton's constant and the cos- 
mological constant in the infrared limit and thus corre- 
spond to those in GR. The above two couplings approach 
non-vanishing constant values in the p <C Mp limit and 
therefore have the expected AS behavior. 

The corresponding RG improved gravitational and cos- 
mological constants obey the following relations. 
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where ^g and ^a are the model parameters determined 
by RG flow coefficients through 



cr; = and t a = —^ 



(9) 



When p — T' 0, they approach the classical values deter- 
mined by observations and thus GR is recovered in the IR 
limit. Conversely, in the extreme UV regime the value 
of G approaches zero, which implies a weakly coupled 
gravitational system at extremely high energy scale. In 
between, we expect a period of sufficiently slow variation 
of A and thus the occurrence of a inflationary phase at 
early times of cosmological evolution. We note that if £^q 
is chosen to be much smaller than unity, one can Tay- 
lor expand the last term of Eq. pi and the simplified 
expression 



A ~ Air + t^U^gGnP 



(10) 



p which can be a function of space-time, varying the La- 
grangian with respect to the metric yields the generalized 
Einstein equation 
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(12) 
where we have introduced the covariant derivative V^ 
and the operator D = — g'^'^V^Vi/. Additionally, consis- 
tency of the Bianchi identity requires the running of the 
cutoff scale to obey the constraint 
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The continuity equation of energy density determines the 
dynamics of matter components and allows derivation of 
the equations of motion by varying the Lagrangian with 
respect to matter fields. The dynamics of this cosmolog- 
ical system are completely determined. 

Inserting the forms of RG modified gravitational con- 
stant G (Eq.(|7|)), and cosmological constant A (Eq.(lO)), 
into Eq.( |13[ ), one can identify the relation between the 
Ricci scalar and the cutoff scale, 
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The original theory may be reformulated as an effective 
f{R) model, where £as is replaced by 
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The subscript ,p denotes the derivative with respect to 
p. In general, the correspondence between the AS gravity 
and f{R) theory holds if the Einstein-Hilbert truncation 
is applied. However, the detailed expression of f{R) de- 
pends strongly on the specific forms of RG functions as 
well as the identification between the cutoff scale and 
Ricci scalar. We refer to [JU] for a general discussion on 
this issue. 



even if p is of order Mp. We will find that this condition is 
necessary in order to achieve a viable inflationary phase. 



III. W INFLATIONARY COSMOLOGY 



B. The /(i?) correspondence 

We consider minimal coupling between the AS gravity 
and the matter field, 



Sas = / d'^x^/^ [£as + -Cm] , 



(11) 



where £m is the Lagrangian of the matter field. As 
the gravitational constant varies along the cutoff scale 



In a realistic cosmological model, the value of Am is 
determined by observations pertaining to late-time ac- 
celeration which are typically of the order O(10~^^^)Afp. 
Therefore, its contribution to the early universe back- 
ground dynamics is totally negligible. For the time being, 
we neglect it and our model reduces to a R^ inflationary 
cosmology J4j. 

In addition, for the matter fleld Lagrangian we focus 
on the SM Higgs scalar. We use the unitary gauge for 
the Higgs boson and temporarily neglect all gauge inter- 
actions. As a consequence, the Lagrangian of the matter 



field is given by 



Cm^-^d^hd^h-V{h)-V,r.t. 



(16) 



where V{h) is the potential of the Higgs boson and Vint 
represents the interactions between the Higgs and other 
particles in the standard model of particle physics. With- 
out considering interactions with other particles, the 
form of the potential is approximately, 



V{h)^\{h^-v^f 



(17) 



in which v is the vacuum expectation value of the Higgs 
boson, with associated Higgs mass, mu = \/2Aw. 



where H = a/ a and the dot denotes the time derivative in 
the Einstein frame. The energy density and the pressure 
in Einstein frame are: 
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By varying the Lagrangian with respect to and /i, we 
obtain the equations of motion for the scalar fields: 
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A. Background equations of motion 



Slow roll inflation 



Now we turn our attention to inflationary solutions. It 
is convenient to perform a Weyl rescaling, 
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where fi is the conformal factor defined by a newly in- 
troduced scalar field as follows, 
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(19) 



As a result, the original AS system is equivalent to 
a two-scalar-field system minimally coupled to Einstein 
gravity without RG running, for which the effective La- 
grangian is given by 

C:,J^-m-SS.-Vi*M). (20) 
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where the potential of the new scalar field takes the form: 

2 
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This potential is sufficiently flat in the regime where ^ 
Mp and migrates into the quadratic form around = 0. 
This scalar field can thus play the role of the inflaton 
under a careful selection of RG running coefficients. 

Substituting the flat Friedmann-Robertson-Walker 
(FRW) metric, ds"^ — ~dt^ + a?{t)dx'^, the Friedmann 
equation becomes 
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Eqs.( 26 ) and ( 27 1 reveal that the inflaton and the Higgs 



fields are coupled and the system is rather intricate. For- 
tunately, the coupling terms are dramatically suppressed 
during inflation due to the slow roll condition. We now 
introduce the series of slow roll parameters 
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(29) 



for a cosmological system of coupled double flelds. We 
note that the potential of becomes very flat when is 
larger than Mp and in comparison the parameters e^ and 
77^ are relatively small. Simultaneously, other parameters 
are also very small due to the suppression by the large 
value of the conformal factor J7. 

As a consequence, under the slow-roll approximation 
the background dynamics are determined by the follow- 
ing solutions: 
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which implies a quasi-exponential expansion at early 
times. Since inflation ends when e^ = 1, the substitu- 
tion of the background solution for 
slow-roll parameter t 
the end of inflation: 



in Eq. ( 28 1 yields the value of 



inEq.(|30|) into the 
at 
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As the contribution of the Higgs field during inflation is 
negligible, the number of e-folding of inflation mainly de- 
pends on through the relation A/" = — /^ Ud(j)/MpU^^, 
which is given by 
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It is easy to achieve Af — 60 ii initially the inflaton is 



placed at 



5.46Mp. Applying the slow- roll condition, 



one obtains the Hubble rate 
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during inflation. 

Moreover, the slow-roll parameters for the Higgs field 
h are automatically small due to the suppression of the 
large value of the conformal factor SI ^ 0(10). Thus dur- 
ing inflation h also varies slowly. As is well known, the 
inflationary phase is an attractor solution in an expand- 
ing universe, and thus it is expected that other matter 
flelds would be dominant in the pre-inflationary phase. 
Specifically, in our model the universe was dominated by 
the Higgs field in pre-inflationary era and at that moment 
the slow roll condition was not satisfied. As a result, the 
parameter r]h can be larger or of order of unity. There- 
fore, we can make use of the relation 77^ ~ 1 to estimate 
the amplitude of the Higgs field at the initial moment of 
inflation, which requires, ry/j = Vj^h/'^H^^'^ ^ 1- There- 
fore, one can estimate the initial amplitude of the Higgs 
at the beginning of inflation as 






Hi 
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At the end of inflation there is no more suppression on 
slow-roll parameters of h and their values become of the 
order of unity simultaneously. As a consequence, the 
amplitude of the Higgs fleld at the end of inflation is 
estimated as 



Vx 



Hr 



(35) 



Combing Eqs.(34| & (351, we easily verify that during 



inflation the Higgs boson is required to satisfy the in- 
equality: 



hf < hi < hi 



(36) 



Note that the Higgs fleld during inflation is not re- 
quired to be less than the Hubble rate since its back- 
ground energy density contributed is dramatically sup- 
pressed by the large value of the conformal factor. Our 
model therefore evades the theoretical constraint sug- 
gested in ^5. We will see that this is key to realizing 
Higgs-modulated reheating. The following numerics ver- 
ify this result explicitly. 

In the above we presented analytic solutions to infla- 
tionary dynamics. We now verify the results with nu- 
merical computations. The results are shown in FigsjT] 
mandS 

FigJT] shows that the Hubble parameter varies very 
slowly in the middle region, which corresponds to the in- 
flationary period. In contrast, the Higgs boson oscillates 
dramatically at the beginning of the evolution, which im- 
plies that the universe is dominated by the Higgs fleld in 




FIG. 1: Evolution of the Hubble parameter H and two scalar 
fields (j) and h as functions of the e-folding number A/". In 
the solutions, the model parameters are, ^g ~ 0.72 and ^a = 
"I^Q-io^^ The parameters of the potential for the Higgs are 
taken as, A — 0.13 and v = 246GeV according to particle 
physics observations. Initial field values are taken as 0^ = 



5.46Mp and ft, 
figure. 



10 Mp. Planck units are adopted in the 




FIG. 2: Evolution of the slow-roll parameters e, e^ and eh as 
functions of the e-folding number A/". The model parameters 
and initial conditions for this plot are the same as those for 
Fig[ll 



the pre-inflation phase. The transition from the Higgs 
dominated pre-inflation phase to the inflation phase fol- 
lows an attractor solution that does not strongly depend 
on the choice of initial conditions. However, this result 
also implies that such a scenario has to meet the big bang 
singularity if one traces backwards in the cosmic evolu- 
tion. From Fig{l] one can see that inflation ends when 
the value of decreases below Mp. The corresponding e- 
folding number is roughly 65. Subsequently, the inflaton 
fleld oscillates around = which corresponds to an IR 




FIG. 3: Evolution of the slow roll parameters r/^, e^ and e^^ as 
functions of the e-folding number A/". The model parameters 
and initial conditions for this plot are the same as those for 
FigB 



fixed point of AS gravity. Therefore, GR is automatically 
recovered at the end of inflation. 

Figj2] and Fig[3] show the evolutions of slow-roll pa- 
rameters defined in Eq.(28) and Eq.(29) along with the 
cosmic expansion. In Figl2l the background slow-roll pa- 
rameter e almost coincides with that for inflaton, e^, since 
inflation is driven by the effects of RG-modifled gravita- 
tional and cosmological constants. Among those param- 
eters associated with the Higgs boson, the value of rjh is 
the flrst to break the slow-roll condition after inflation. 
Consequently, the method of determining the value of h 
during inflation by requiring e^ ~ 1 is reliable. By sub- 
stituting the parameter choices provided in Figjljinto the 
expression Eq.(36), we flnd hj ~ 10~^, which is consis- 



tent with the numerical result shown in the lower panel 
ofFigjl} 



C. Higgs dependent decay after inflation 

After the inflaton field rolls below the critical value 
(j)f , it starts to oscillate around the IR point which corre- 
sponds to the GR limit. One can thus perform a Taylor 
expansion of the potential around (f> = 0: 
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up to the (jP' order. In the above expression, we have 
introduced an effective mass for the inflaton defined by 
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(38) 



The last term of Eq. ( 37 1 shows that cj) interacts with other 



particles through the expansion of the conformal factor. 



Thus, if Vint contains interactions of the Higgs boson with 
other particles, the same interactions provide channels for 
the inflaton to decay into them and the corresponding 
decay rate is expected to be a function of the Higgs field 
value. Before the inflaton decays, the evolution of the 
universe is dominated by the mass term and thus the 
background equation of state is effectively Wm = 0. 



In the current case, the last term of (37) is respon- 
sible for the Higgs-dependent inflaton decay. Follow- 
ing Ref.[53], one can generally take the following Higgs- 
dependent interactions: 

VL ^ Va[h)c^Ma + Ma{h)cj)xl + 9a{h)4>''xl , (39) 

where Xa and f/'a are the scalar and spinor fields which 
constitute radiation in the early universe; the subscript a 
represents the species of particles. To achieve the mod- 
ulated reheating scenario in our case, the coupling con- 
stants Ha, Ma, and ga must be functions of the Higgs 
field. Under this assumption the decay rate of the in- 
fiaton to the lowest order in coupling constants is given 
by 
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(40) 



where the quadratic potential for </> has been applied. 

Having obtained the decay rate of the inflaton, we now 
calculate the time when the inflaton fleld decays com- 
pletely, and the universe enters the radiation dominant 
phase. In the following, we examine a specific example 
to illustrate the possibility of this decay process. We 
consider 
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which could arise from the term Vint't''^ / M^ appearing in 
the last term of Eq. ( 37 ) . The corresponding decay rate 
is given by 
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At the moment of the phase transition from inflaton dom- 
ination to radiation domination, we have the condition 
i7 ~ r. As we will show, this is the spatial hyper-surface 
where the primordial curvature perturbation can be cal- 
culated. At the reheating surface, p^ ~ SMpV'^. By 
making use of Eq. (351, we obtain the value of inflaton 
decay. 
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(43) 



at the reheating surface. 



IV. MODULATED REHEATING VIA HIGGS 
BOSON 

In this section, we briefly review the standard calcula- 
tions of the primordial power spectrum, the bispectrum 



and the trispectrum for the mechanism of modulated re- 
heating, with the assumption that the inflaton decays on 
a spatial hypersurface with a variable local decay rate. 
Afterwards, we will apply this mechanism to our model 
in which the decay rate is a function of the Higgs boson 
and then study its cosmological implications. 



A. Field fluctuations during inflation 

We analyze primordial perturbations in a double-field 
inflation model involving kinetic couplings. We refer to 
[5H - I55] for earlier studies of inflation models in terms 
of kinetically mixed double fields and [571 ISH] for the 
paradigm of double-field inflation. Also, the topic of pri- 
mordial perturbations in multiple-field inflation models 
was recently reviewed in [SSI EO] ■ 

During inflation, it is convenient to decompose the field 
space of our model to directions parallel and orthogonal 
to the trajectory of background evolution. Along these 
two directions, one can define the adiabatic field, a, and 
the entropy field, s, as follows, 
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s ~ — sin 6 



n ^ cos 9h , 



where the rotation angle is given by 
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After that, we perturb the metric and fields up to linear 
order. One can introduce the field fluctuations along the 
adiabatic and entropy directions as follows. 
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Neglecting interactions between these two modes, one 
can solve for the amplitude of entropy perturbation, 
Ss^ = 7f*/27r, during inflation and therefore the field 
fluctuation for the Higgs boson, which is given by 



6h^ — fl^Ss^ — Q, 



2^ ' 



(49) 



where the subscript * denotes the moment of Hubble 
crossing. 

B. Higgs-modulated reheating 

The generation of curvature perturbation via Higgs- 
modulated reheating in a canonical model was recently 
suggested in j48l 09] , and soon it was pointed out in [47] 
that the Higgs-dependent interaction potential of the in- 
flaton would be severely constrained by an upper bound 
on the value of the Higgs during inflation. 



In the present paper, we extend the paradigm into the 
non-canonical model under consider. We show that this 
upper bound can be greatly relaxed by the relatively large 
value of the conformal factor and thus the corresponding 
parameter space is dramatically enlarged. Our scenario 
is easily extended to non-minimal inflation models. 

In the treatment of local non-Gaussianity, the curva- 
ture perturbation can be expanded order by order as fol- 
lows, 



C(x) = a{x) + ^hUU^) + ^5nlCi'(x) + 0{Ct) 



E 



C,n{x) 



(50) 



where Ci is the Gaussian fluctuation and Cn are the non- 
Gaussian components of order C". The relation between 
C„ and the non-Gaussian parameters yields the following 
non-Gaussian estimators, 



6Ci 



25 C3 



(51) 



The correlation functions are defined as 

2 

(C(ki)C(k2)) = (2^)3p(fcl)^'(Ekn), 

n=l 

3 
(C(ki)C(k2)C(k3)) - (2^)3i3(ki,k2,k3)(53(^k„) , 

n=l 

4 

(C(ki)C(k2)C(k3)C(k4)) = (27r)3r(ki,k2,k3,k4)<53(5]k„) 

n=l 

(52) 

where -P(fci) is related to the dimensionless power spec- 
trum as 



Mki 



fc3 



-P{k,) 



(53) 



Inserting the ansatz Eq.(50l into Eq.(52), one can relate 



the bispectrum B and the trispectrum T with P as fol- 
lows, 

B(ki, k2, kg) = ^/nl [P{ki)P{k2) + 2 perm.] , (54) 
54 



r(ki,k2,k3,k4) = -ffNL [P{k{)P{k2)P{h 

25 



3 perm.] 
TNL [P(fci)P(fc2)-P(|ki + kgl) + II perm.] . (55) 



Note that if we neglect the nonlinear perturbations in- 
duced by self-interactions during inflation, which will be 
treated in the subsection |IVD[ there exists, in this case, 
a simple relation tnl — (36/25)/^l. 

In the modulated reheating scenario, the decay of the 
inflaton occurs on a spatial hyper-surface with a varying 
local decay rate F, which is assumed to be a function of 
the Higgs boson in our model. Thus, the local Hubble 



parameter on the slice of modulated decay satisfies the 
condition H = T{hy. On super-Hubble scales, the cur- 
vature perturbation arisen from modulated decay can be 
written as 



Ux) 



-Gi 



Sh 
~h 



1 



rOo 



763 



,(56) 



where the subscript d denotes the moment of modulated 
decay. In our model the </> potential is dominated by 
its mass term after inflation, as shown in Eq.(37). We 
therefore obtain the coefficients as follows, 



61 = 


hV.h 

6 r ' 


62 = 


" 6 ^ r 


63 = 


' 6 ^ r 



-), 



r,,r 



,/iJ- ,hh 
"f2~ 



p3 



(57) 
(58) 
(59) 



At linear order in the curvature perturbation, the co- 
efficient Gi is typically of the order 0(1) and thus C is 
mainly determined by 5h/h at the moment of modulated 
decay. In the conventional scenario of modulated reheat- 
ing, one can approximately take 5h£, to be the amplitude 
of entropy field at the moment of Hubble-crossing dur- 
ing infiation. However for the model of AS inflation, the 
Higgs boson and the inflaton are coupled through a con- 
formal factor in front of its kinetic term. Making use of 



Eq.(49l, one finds. 



5h 



D 






(60) 



Moreover, the value of the Higgs at the slice of modu- 
lated decay can be related to the Hubble-crossing value 
by introducing a general function 



^D 



g{K) , 



(61) 



where its detailed form is determined by the explicit po- 
tential of the entropy field. For example, in the model 
under consider, g{h^,) ex /i* ^ hj. For simplicity, we 
assume that gih^,) is linear. 

As a consequence, the power spectrum of curvature 
perturbation due to modulated reheating is given by 



Pc 



a 



"■D 



eln^ 



_Hl_ 

^47r2;i2 



(62) 



where we have applied the field fluctuation Eq.(60) to 



obtain the second, approximate equality, expression. 



^ Within the framework of the multi-field inflationary cosmology, 
there exist many interesting scenarios for generating primordial 
curvature perturbation based on different choices of decay slices, 
such as the modulated curvaton decay mechanism I61H63I and 
the uniform curvaton decay mechanism [64]. All these scenarios 
are well established based on the validity of the SN formalism. 



C. Observables at linear order 

If we further neglect the variation of the Hubble param- 
eter during inflation, then we can obtain an approximate 
power spectrum from Higgs-modulated reheating. 



P, 



A 



Ch 



02 



(63) 



by inserting the approximate relation Eq. ( 34 1 . Note that 



the usual decay rate is a power-law function of the Higgs 
boson such as that in Eq.(42), considered in the previous 
section. Thus Qf is typically of the order 0(0.01 ~ 1). 
The coefficient flo is totally determined by (/i^i as pro- 
vided by Eq.(|43]), and it therefore depends only on A 
and k; numerically flD is of the order 0(1). Finally, we 
note that the power spectrum generated from the Higgs- 
modulated reheating is determined by the Higgs coupling 
A, the interaction coupling k, and the conformal factor 
during inflation 0/ (or, cquivalently, the value of inflaton 

In addition to the curvature perturbation generated by 
modulated reheating, there exists the intrinsic curvature 
perturbation due to the inflaton fluctuation, which takes 
the form: 



Pc- 



Hf 



87r2eM2 



(64) 



It is convenient to define a Higgs-to-curvature ratio 



P 



qh 



Xh 



exeln^M 



Pa+Pc, exelnlM^ + n^jHf 



(65) 



to characterize the relative contribution of Higgs fluctua- 
tions. If the main contribution to generating primordial 
curvature perturbation is due to the modulated reheat- 
ing, then we expect qh — I. By choosing a group of 
values for the model parameters such as that provided 
in FigII|and the decay rate given by Eq.(42), one finds: 
e ~ 10-*, ni - 10, Hd - 1, e? - 0.1 and A - 0.13. 
Under this particular parameter choice, we find that the 
mechanism of Higgs modulated reheating dominates as 
long as Hj < 10~^Mp without any fine-tuning. 

These two power spectra actually show different signa- 
tures on their spectral indices. Specifically, their spectral 
indices are given by 



n^_^ - 1 = -6e + 27^0 = -6e 



2C/.. 



\^ 



1 



-2e + 2r]h = -2e 



3i/2 ' 



(66) 
(67) 



which are calculated at the moment of Hubble-crossing. 
In addition, the primordial tensor perturbations are only 
dependent on the inflationary Hubble parameter, whose 
spectrum is given by 



Pt 



2IP_ 

7r2M2 ' 



(68) 
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As usual, the spectral tilt is given by 
riT = — 2e . 



(69) 



In the modulated reheating scenario the conventional 
tensor-to-scalar ratio rx is now defined as 



Pt 



tt 



Pc.+Pc, 



16e(l - qh) , 



(70) 



which indicates that the amplitude of primordial gravita- 
tional wave is doubly suppressed in the Higgs-modulated 
reheating mechanism since both e and 1 — q^ are small 
quantities. 



cess of post-inflation modulated reheating. In this era 
the Higgs field is considered Gaussian while the non- 
Gaussianity is induced by the nonlinear conversion from 
dh to ^. One can insert the second and third order cur- 



vature perturbations in Eq.(56) into the non-Gaussian 



estimator Eq.(51) and obtain this part of the "universal" 



nonlinearity parameters: 



j-local 
/NL,un 



^local 
yNL,un 




(71) 
(72) 



D. Non-Gaussianities 

In contrast to the prediction of a canonical single-field 
inflation model ^SF, a salient feature of the modulated 
reheating mechanism is that sizable amplitudes of pri- 
mordial non-Gaussianities can be obtained under suit- 
able parameter choices. In this subsection we study the 
curvature perturbation beyond linear level. As we ex- 
pected, the curvature perturbation is mainly sourced by 
the Higgs fluctuations. For the time being, we ignore 
the nonlinear effects of inflaton, which are generally sup- 
pressed by slow-roll parameters. For the nonlinear fluctu- 
ations seeded by the Higgs fluctuations, there exist two 
categories of seeds, with one being proportional to the 
connected correlators of the Higgs and the other being 
an intrinsically non-Gaussian distribution |49j . 



1. Non-Gaussianities from modulated reheating 

The first type of non-Gaussianities originates from the 
field fluctuations at super-Hubble scales during the pro- 



which are of local type. 

In particular, for the interaction term considered in 
(41 1, the decay rate is proportional to h^ and q/,, ~ 1 can 



be obtained under a reasonable set of values of model 
parameters. As a consequence, one obtains /nlIIii — 5/2 
and ^nl'uii — 25/3. These nonlinear parameters are 
sizable when compared with those in slow-roll inflation 
models, but the corresponding non-Gaussianities are still 
difficult to test observationally. 
2. Non-Gaussianities from Htggs Self -interaction during 
inflation 



The second type of non-Gaussianities originates from 
the non-quadratic potential of the lighter field, which in 
our model corresponds to the Higgs potential: V{h) ~ 
A/i'*/4. In fact, this self-interaction of the scalar field 
can also generate primordial non-Gaussian fluctuations 
during inflation. 

Following [49] (see also [66]), the n-point correlation 
function of 5h is evaluated by 



('5/iki(T),5/ik.(T) • • • 5h^Ar)) = -«(l / adT'[6h^,{T)5h^,{T) ■ ■ ■ dh^^r), H^!;il{h(T') 



(73) 



where H^^t ^^ ^he n-th order interaction Hamiltonian. 
Here by n-th order we mean the part of Hint that is of 
the order 0{Sh^). In our model the Higgs field is con- 
formally coupled to the infiaton due to the RG running 
gravitational constant. The corresponding field fiuctua- 
tion is expressed as 



Shv 



inn 



(1 4- ikT)e 



-ikr 



(74) 



during infiation. In addition, the interaction Hamiltonian 



of the Higgs field takes the form 

d3a;a3O-''(0)^y(")(/i)<5/i" , (75) 



where V'^^\h) = d^^V/dh^ is the n-th derivative of the 
potential V{h) with respect to the field h, and V{h) as 
well as ri((/)) are given by Eq.(17) and Eq.(19), respec- 
tively. 

We perform the integrals appearing in the correlation 



10 



functions and find 



((5ft,ki(5ft.k2 ■ ■ -Shk 



{n,H, 



\2n-4 



VJ^'^K^ 



nLi(2fcf) 

x(2^)3/„(ki,k2r 



s'ij:^. 



4=1 



bispectrum) and ki <^ fc2,fc3,fc4 (for trispectrum). The 
same scenario in the framework of GR was discussed in 
|49) . Here we directly calculate the correlation function 
of the curvature perturbation and then derive the non- 
linearity parameters: 



(76) 



where a kernel integral function has been introduced as 
follows, 



2Re 



dr' 



W{1 - ikiT)e 



iKr 



i=l 



(77) 



with K defined as i^ = X^ILi ^i- 

For 3- and 4-point correlation functions which are of 
observable interests, we identify, according to the defini- 
tions in Eq.(52|, the following expressions: 



i?,V""(ki,k2,k3) 



nli(2fcf) ' 
4nr=ifcf 

2 



2 l^i<j '^i'^j 



9 



K-^ 



(78) 



and 



T> 



Sh 



(ki,k2,k3,k4 



-f2n^- 



2E 






■{ie + Nk) 



E.fcf 



K3 



where je — 0.58 is the Euler-Masheroni constant and 
Nk is the e-folding number for the perturbation mode 
with a fixed K crossing the Hubble radius until the end 
of inflation Tend- As introduced in previous section, the 
subscript "*" indicates the values at Hubble crossing. 

We first calculate the non-Gaussianities of equilateral 
type. One can estimate the correlation functions Bgf^""' 
and Tgf^^" under the particular limit that all the fc^'s are 
of the same value. As a result, substituting Eqs. (78 1, 



(79) into the expressions Eq.(52) yields the nonlinearity 



parameters of equilateral type as follows. 



/.cquil 
•/NL.int 



equil 

^NL.int 



qt[NK + iE-i 



25Xhl 3/ ^ 169\ 



(80) 
(81) 



Next we study the non-Gaussianities originated from 
the self-interaction of the Higgs field during inflation in 
the squeezed limit where we assume ki <C ^2,^3 (for 



j-local 
/NL,int 



local 
ffNL.int 



5Xhl 



2bXhl 



Nk +1e 



qt NVk + 7E - 3 



(82) 
(83) 



(79) 



From the above results, we can immediately see 
that the primordial non-Gaussianities due to the self- 
interaction of the Higgs field during inflation are negative. 
This is a novel feature in the Higgs modulated reheating 
scenario. A similar feature was observed in |49] in the 
framework of standard GR, but in our model the ampli- 
tude of nonlinearity parameters involve a new parameter 
which is the conformal factor J7. 

Consider for the moment the primordial curvature per- 
turbation due solely to the modulated reheating, g^ ~ 1. 
By choosing a group of reasonable values for the model 
parameters such as that provided in the previous section, 
we find e - IQ-*, ^i - 10, 6? ~ 0.1 and A ~ 0.13. In ad- 
dition, there is a theoretical lower bound: hi > 2Hi/y/X. 
By assuming Nk ^ 50, one obtains /NL^int ^ ~10. We 
see that this particular parameter choice appears to be 
incompatible with the newly released Planck data. This 
points to the necessity of performing an analysis of the 
observational constraints on our model. This is the main 
content of the next section. 



V. CONSTRAINT ON MODEL PARAMETERS 
BY PLANCK 

Recently the Planck mission has released data on CMB 
anisotropy. The results highly constrain cosmological 
parameters with unprecedented accuracy. Specifically, 
the amplitude and spectral index of primordial curvature 
perturbation are determined to be lO^P^; = 2.23 ± 0.16, 
n, = 0.9603 ± 0.0073 (68% C.L.) at the pivot scale 
k = 0.002Mpc~^ [S]. Moreover, there is no significant 
evidence for primordial curvature perturbation deviating 
from Gaussian distribution. In particular, the bounds on 
nonlinearity parameters are quoted as: /^l^' = 2.7±5.8, 
^oquai ^ _42±75 (68% C.L.) [1]. In addition, the upper 
bound on the tensor-to-scalar ratio is given by r^ < 0.11 
at two-sigma level. 

In our model of RG improved Higgs modulated reheat- 
ing, there are eight model parameters. Among these pa- 
rameters, Hif, flif and e are associated with the back- 
ground model; /i* and A are related to the details of the 
Higgs model; and 9i, 02 and 83 are determined by spe- 
cific forms of the decay process, respectively. To be more 
explicit, A is basically constrained by particle physics ex- 
periments such as those at LHC, which have determined 
that A ~ 0.13. At present the Planck data has not yet 
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imposed strong constraints on the tri-spectrum and thus 
03 is free. We therefore only need to analyze the com- 
bined constraints on the remaining parameters. 
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FIG. 4: Observational and theoretical constraints on /i* and 
Hf The viable parameter space is within the red region. As 
Qh increases from zero to unity, the color reduces to a lighter 
tone, which corresponds to a more severe fine-tuning problem 
in this limit. Planck units are adopted in the figure. 

Allowing all plausible values for e and Oi, we obtain a 
combined constraint on the inflationary Hubble rate iJ* 
and the amplitude of the Higgs field ft.* as depicted in 
Fig|4] One can read from the figure that H^, is tightly 
constrained to be of the order O(10~'''^Mp) and h^, is al- 
lowed to vary between 10"'* Mp and 10~^Mp during infla- 
tion. In Fig|4]the color of the region reduces to a lighter 
tone as qh increases from zero to unity, which corresponds 
to a more severe flne-tuning problem towards the qh = 1 
limit. This indicates that the case for the pure Higgs 
modulated reheating, where q^ — 1, is already ruled out 
by the Planck data. The combined scenario of Higgs 
modulated reheating and inflationary perturbation, how- 
ever, survives. 



VI. CONCLUSION 

In this paper, we studied the RG improved inflationary 
cosmology where the decay rate of inflaton is modulated 
by a second scalar field, which we identified as a Higgs. 
In this model the background evolution is driven by a 
RG running cosmological constant and gravitational con- 
stant with their RG flows satisfying the AS behavior. By 
choosing Einstein-Hilbert truncation, we flnd this model 



is classically equivalent to a model of R^ gravity. The 
elegant property of this model is that it can give rise to 
a sufficiently long inflationary phase at high energy scale 
and smoothly exit to standard GR after inflation. More- 
over, an RG running gravitational constant can assist a 
second scalar field to vary slowly without an extremely 
fiat potential since the slow-roll parameters associated 
with this field are greatly suppressed by a large value of 
the conformal factor. As a consequence, this scalar field 
seeds isocurvature perturbations during inflation which 
can be converted into primordial curvature perturbation 
under a suitable mechanism. We consider this mecha- 
nism as the process of modulated reheating. 

Based on this mechanism, we performed a detailed 
analysis on the power spectrum and non-Gaussianities 
of primordial cosmological perturbations. We then con- 
fronted our model with the recently released Planck data 
and concluded that a viable parameter space exists, al- 
though it is highly constrained. In order to be consistent 
with observations, we are forced to combine the curvature 
perturbations generated by the Higgs modulated reheat- 
ing with the fluctuations seeded by the inflaton itself. 
Although this model suffers from the fine-tuning prob- 
lem, the scenario under present study points to a new 
possible connection between particle physics and early 
universe cosmology. 

We conclude by mentioning that the mechanism of 
Higgs modulated reheating can be generalized to an arbi- 
trary non-minimal inflationary model or a model of f{R) 
inflation. By relaxing theoretical requirements of AS 
gravity, the parameter space available to such a mech- 
anism is increased. This topic is of phcnomenological 
interest. 
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